Serum concentration of polychlorinated biphenyls and the risk of type 2 diabetes: a 10-year follow-up historical cohort study

This study investigated the association between serum concentrations of Polychlorinated Biphenyls (PCBs) and the risk of type 2 diabetes within the general population. A ten-year follow-up historical cohort study was conducted during 2009–2019 as part of the Bushehr MONICA cohort study in Iran. Of 893 non-diabetes participants at base line, 181 individuals were included in the study. The concentration of nine PCB congeners was measured in individuals’ serum samples at baseline, and the risk of type 2 diabetes was determined based on fasting blood sugar at the end of follow-up. Multiple logistic regression models were used to assess the study outcomes after adjusting for covariates. This study included 59 diabetes individuals (32.6%; mean [SD] age: 58.64 [8.05]) and 122 non-diabetes individuals (67.4%; mean [SD] age: 52.75 [8.68]). Multivariable analysis revealed that a one-tertile increase (increasing from 33rd centile to 67th centile) in Σ non-dioxin-like-PCBs (OR 2.749, 95% CI 1.066–7.089), Σ dioxin-like-PCBs (OR 4.842, 95% CI 1.911–12.269), and Σ PCBs (OR 2.887, 95% CI 1.120–7.441) significantly associated with an increased risk of type 2 diabetes. The strongest association was obtained for dioxin-like PCBs. The results highlight a significant correlation between PCB exposure and an increased risk of type 2 diabetes. The evidence suggests that additional epidemiological studies are necessary to clarify the link between PCBs and diabetes.


Study population
We conducted a historical cohort study utilized baseline data from phase II of the PIB-MONICA study.Participants who were in phase II of the study between 2009 and 2010 and followed up for 10 years, were invited to this study conducted between 2019 and 2020.All non-diabetic individuals, based on their registered data at the end of phase II, were identified (n = 893 out of 1092).We were informed about their condition and current use of insulin or anti-diabetic drugs through a telephone interview.All non-diabetic individuals residing in Bushehr, Iran, were eligible for this study (n = 515), out of which 70 reported having diabetes and 445 reported no diabetes based on the phone self-declaration.The remaining participants (n = 378) were excluded for various reasons, such as death, relocation, non-response due to a phone failure or changed numbers, or unwillingness to participate in the study.The number of samples was limited to the number of eligible people from the baseline who were willing to participate in this study.For calculating the study sample size, we used Codru et al. 21, a cross sectional PCBs evaluation in diabetic adult Native-Americans, as a reference.We calculated the sample size with a power of 80%, alpha = 0.05, the odds ratio (OR) of 3.9, and the prevalence of exposure with diabetes in case group of 0.33.The planned sample size for case group was 53 and considering r = 2 (the ratio of control to case), the sample size of the control group was also calculated.Taking into consideration the 15% drop-out rate, we indicating that a total of 183 participants would be required.Since the sampling period coincided with the COVID-19 lockdown, it affected the participation rates, especially in the case group, and we could not achieve the desired number of cases.Ultimately, 181 individuals (59 people with diabetes and 122 non-diabetes) were included in the study.
Subsequently, blood samples were taken from all participants to measure some biochemical characteristics, including fasting blood sugar (FBS), total cholesterol (TC), high-density lipoprotein (HDL) cholesterol, and triglycerides (TG).Low-density lipoprotein (LDL) cholesterol levels were calculated using the Friedewald formula 27 .Additionally, other anthropometric indices such as blood pressure (measured twice after 15 min of rest in the right arm), height, weight, and waist and hip circumference were measured and recorded using a modified retrospective data collection form, which was based on the standardized WHO MONICA Project questionnaire 25,26 .

Outcome assessment
The primary outcome, the risk of type 2 diabetes, was assessed through fasting blood sugar (FBS) serum levels.For this purpose, we considered the following definitive criteria for the participant with diabetes: having fasting blood sugar (FBS) ≥ 126 mg/dl [28][29][30] or consuming anti-diabetes drugs.

Covariates
Data were collected via face-to-face interviews with a professional interviewer using the data collection form that included demographic information, tobacco use, physical activity, hypertension, dyslipidemia, body mass index (BMI), other diseases status, occupational exposure status, presence in cities with air pollution, aquatic food consumption pattern, and supplements status such as vitamin D and omega3.Physical activity scores were calculated using the short version of the International Physical Activity Questionnaire (IPAQ).Subsequently, the scores were classified into three levels (high [HEAP active], moderate [minimally active], and low activity [inactive]) 31 .According to NCEP ATP-III hyperlipidemia criteria 32 , the participants with dyslipidemia were identified.In addition, participants with hypertension were identified based on WHO criteria 33 or on a history of consuming antihypertensive drugs.

Sample collection
The PCBs serum concentration was measured to describe the environmental exposure in this cohort study.To evaluate the exposure level of 181 participants to PCBs, their blood serum samples, stored at − 80 °C for the past 10 years, were utilized.These serum samples were collected and preserved at the end of phase II of the PIB-MONICA-S in 2009.The concentration of nine PCB congeners was determined based on the International Union of Pure and Applied Chemistry (IUPAC) nomination, including 28, 52, 99, 101, 114, 118, 138, 153, and  180.Then, the PCBs exposure was categorized into three clusters of total PCB concentration (ΣPCB: the sum of the nine studied PCB congeners), NDL-PCBs concentration (the sum of seven PCB congeners: 28, 52, 99, 101, 138, 153, and 180), and DL-PCBs concentration (the sum of two PCB congeners: 114, and 118).

Serum extraction
The laboratory extraction of PCBs was conducted in a certified private laboratory accredited by the Tehran University of Medical Sciences, employing solid-phase extraction (SPE) as the extraction technique.The laboratory extraction of PCBs was conducted in a certified private laboratory accredited by the Tehran University of Medical Sciences, employing solid phase extraction (SPE) as the extraction technique.The sample concentration and analytical procedure were guided by the methodology described in the study by Lin et al. 34 .
Solid-phase extraction, a modern extraction technique utilized in this study, encompasses four steps: conditioning (to prepare the stationary phase), sample application (to add the sample into a column containing the stationary phase), washing (to wash the column with a suitable solvent to remove interfering substances), and finally elution (to wash the column with a suitable solvent for analyte extraction).The procedure employed a vacuum manifold, nitrogen gas, carbon-bonded cartridges, and a glass connector for solvent concentration and drying.Also, Dichloromethane was used for eluting the analytes under a vacuum following SPE conditioning steps.

Gas chromatographic conditions
The PCBs congeners were analyzed and measured using an Agilent 6890-5973 gas chromatography-mass spectrometry (GC-MS) equipped with an HP-5MS column (0.25 µm film thickness, 0.25 mm, i.d., 30 m).The initial temperature of the GC oven column was 90 °C, maintained for one minute, then increased to 150 °C at a rate of 50 °C/min for one minute, and finally increased rapidly to 330 °C at a rate of 8 °C/min, where it was held constant for three minutes.The injector and detector temperatures were adjusted to 250 and 150 °C, respectively.Helium was applied as the carrier gas at a fixed flow rate of 1.0 ml/min, and the transmission line temperature was set at 280 °C.This method enabled the detection of desired analytes in nanograms per liter.
Validation procedure (quality assurance/quality control (QA/QC)) A calibration curve was plotted with concentrations of 0.3, 1.5, 5, 15, 25, and 50 ng/L for the mixture of targeted PCB analytes (standard solutions).Linear regression was applied to plot the calibration curve area against concentration.Pearson correlation coefficient was calculated and obtained from 0.96 to 0.99 for various analytes.Moreover, qualitative identification of analytes was done based on the retention time (RT), ion spectrum, and ion ratio.The filament was turned off during the initial 2 min of detection to prevent solvent interference.Hexachlorobenzene used as the internal standard with a concentration of 0.1 μg/liter, which was added to the samples after dissolving in methanol.To check the accuracy of the quantitative assessment, within-day and between-day precision was evaluated in triplicate over three separate days for low, medium, and high standard PCB analyte concentrations (0.5, 7.5, and 35 ng/l).The limit of detection (LOD) and the limit of quantification (LOQ) were obtained equal to 0.1 and 0.3 ng/L, respectively, for each measured congener.

Identification and quantification of PCBs
PCB congener concentrations were measured based on a wet weight basis.Due to their lipophilic nature and tendency to associate with blood lipids, a concentration on a lipid-adjusted basis was used 35 .Total serum lipid content was calculated by summing the triglycerides (TG) and total cholesterol (TC) according to the Phillips formula 36 , with values expressed as nanograms or picograms per gram lipid (ng/g lipid or pg/g lipid).

Statistical analysis
The normality of the data is determined using the Kolmogorov-Smirnov test and the histogram.Homogeneity of variance was assessed using Levene's Test for Equality of Variances.Data were examined using descriptive analysis, which included the mean (SD), median (interquartile range (IQR)), and frequencies (percentages).The relationship between quantitative variables and the dependent variable was examined using independent t-tests, while chi-square tests were used for categorical variables.Simple logistic regression models tested between type 2 diabetes and demographic, anthropometric, biochemistry characteristics, diet and supplement information of participants, and other diseases separately to identify the potential confounders.Finally, age, waist circumference, family history of diabetes, and dyslipidemia were considered confounding variables.Due to a relatively high loss to follow-up, odds ratio (OR) approach were used instead of relative risk (RR).Multiple logistic regression model was used to evaluate the relationship between PCB serum concentrations with type 2 diabetes.To control the effect of confounders, other covariates (age, waist circumference, family history of diabetes, and dyslipidemia) were included in the model.The obtained data from multiple logistic regression models were demonstrated using OR with confidence interval (95% CI).The Hosmer-Lemeshow test was used to check the goodness-of-fit for the regression model.The model with the highest steps has the best fitness.The level of statistical significance in this analysis was 0.05.In addition, if the reported concentration of the PCB congeners was within the limit of detection or below, half of the detection limit was recorded instead (LOD/2).IBM SPSS Statistics 19 was used for all statistical analysis.

Ethical considerations
The study protocol and all procedures were approved by the Tehran University of Medical Sciences Ethics Committee (# IR.TUMS.SPH.REC.1398.321).All participants signed and confirmed the written informed consent.Moreover, Participation in the study was voluntary.

Baseline characteristics
We studied a total of 181 participants during ten years of follow-up (Fig. 1), of whom 59 (32.60%) were categorized as having type 2 diabetes; most diabetic participants were females (31 of 59 [52.50%]), and 54 of 59 [91.50%] were married.The mean (SD) age of participants was 54.67 (8.9) years, significantly higher in diabetic than non-diabetic participants (p < 0.05).Other characteristics of participants including clinical and laboratory details during follow-up, are summarized in Table 1.Participants at risk for type 2 diabetes were more likely

Chemical concentration in participants' blood
Altogether, 181 serum samples from the beginning of the follow-up period (in 2009) were analyzed to measure the concentration of PCBs.The lipid-adjusted serum concentrations of ΣNdl-PCBs, Σdl-PCBs, and ΣPCBs were significantly higher in the diabetes group (p < 0.05) (Table 2, Fig. 2A).PCB congeners 153 and 138 exhibited the highest concentrations.Also, they were reported with approximately equal proportions between both groups (PCB138: 40% in diabetes vs. 40.50% in non-diabetes and for PCB153: 50.50% in diabetes vs. 52.30% in nondiabetes).Except for PCB 138 and 153, all tested congeners had lower values than the detection limit in some samples.The concentration of total PCBs was detected in the serum of all participants so that the concentration of any sample was not reported as zero.

Regression models
The logistic regression models for PCBs are shown in Table 3 and Fig. 2B.

Discussion
In this study we aimed to assess the relationship between the serum PCB levels in individuals with type 2 diabetes and those without.This study utilized the high-quality MONICA cohort data with a ten-year follow-up period, incorporating individuals without diabetes at baseline.Based on this study result, the detected serum PCB concentrations ranged from 0.37 to 15.68 ng/g lipid, with a median value of 4.08 ng/g lipid (interquartile range [IQR] 2.51-6.15).Also, all serum samples exhibited some level of PCB concentration.When compared with findings from global studies, the results revealed lower concentrations than those reported in Pakistan 37 , China 38 , and Japan 39 .These relatively low levels suggest a lesser exposure of the study population to PCBs; however, these minimal concentrations do not diminish its importance.In this regard, scientists also pointed out that exposure to POPs, even at low concentrations, can have significant human health outcomes, including an elevated risk of diabetes 5,40 .Neglecting the population groups with lower exposure risks, compared to those with occupational or accidental exposure, might mask or overlook potential hazards 40 .On the other hand, the results are within the range of PCB concentration detected on occupational exposure groups in Tehran study 41 .Interestingly, the findings indicated higher concentrations than a study involving pregnant women in Tehran, likely attributable to Bushehr's greater consumption of fish and aquatic foods compared to the national average 42 .Food sources are estimated to account for approximately 90% of PCB exposure portion in Iran, as reported by Ravanipour et al. 43 .Also, PCB153 and PCB138 were detected in all serum samples (with a detection rate of 100%).Similarly, some studies reported these congeners as the www.nature.com/scientificreports/most abundant congeners 15,38 .It is generally accepted that the three constituents of PCBs (153 + 138 + 180), on average, about 61% of the PCB burden in the human body 44 .This study analyzed nine main congeners, categorized into dioxin-like and non-dioxin-like PCBs, and included in the models after adjusting for confounders.The literature displays variability in the number and types of congeners analyzed and their classification into varied subcategories for assessing their association with type 2 diabetes.Some studies have focused on individual congeners (e.g., PCB 153) 17 or a combination thereof, but most aggregate them into the total PCBs (ΣPCBs) metric 5,18 .Accordingly, the number of congeners used to calculate the total PCBs varied widely.The Stockholm Convention recommends measuring six to seven indicator PCBs (PCB-28, 52, 101, 138, 153, 180, and 118) due to their higher concentration in food, the environment, or human tissues/fluids 8 .These noteworthy methodological differences may influence the observed relationship with type 2 diabetes.
A significant positive correlation was identified between total PCBs, non-dioxin-like PCBs, and dioxin-like PCBs with type 2 diabetes; the strongest association was observed with dioxin-like PCBs in this study.While some studies support the findings, others report no significant correlation.Airaksinen et al. (2011) confirm a positive association between type 2 diabetes and exposure to PCB153 in a general urban population of the Helsinki Birth Cohort Study 17 .Philibert et al. 18 also reported a positive association between total PCBs and type 2 diabetes in 101 Canadian participants.In contrast, Jorgensen et al. 16 in Greenland examined the relationship between 13 of the most common PCBs after adjusting for age, sex, waist circumference, hereditary history, smoking, alcohol consumption, and education level, but no correlation was found.Also, Zani et al. 15 found no association between total PCBs and diabetes in their study of an urban area in northern Italy where a PCB chemical plant had previously caused severe contamination.Similarly, in a cohort study by Magliano et al. 2 with nine years of follow-up in France, no association was found between exposure to ΣPCB and diabetes.Although most of these studies have been performed as a case-control or cohort study, the results of meta-analysis studies obtained from the aggregation of such articles have not shown much certainty in their conclusions.These discrepancies underscore the complexity of determining the impact of PCB exposure on diabetes risk.A meta-analysis study conducted by Tang et al. 45 showed a possible association between PCB153 and type 2 diabetes.The need to conduct smaller studies with more inconsistent results is emphasized in future studies in order to minimize the potential impact of publication bias.Conversely, a meta-analysis study conducted by Song et al. 24 emphasize the urgent need to investigate the relationship of endocrine-disrupting chemicals (EDCs) with diabetes risk, mainly through large prospective studies.
Moreover, today, the main mechanism of PCBs association with diabetes is controversial and still unclear.Various in vitro and in vivo studies have addressed this issue and mentioned some possible mechanisms, including mitochondrial dysfunction, endocrine disruption (through depletion of pancreatic β cells), and the effect on glucose metabolism 46 .In addition, the various toxic effects of PCBs are related to the chlorine site on the phenyl ring.The high ability of some PCB congeners for binding to aryl hydrocarbon receptor (AhR) causes toxicological properties similar to polychlorinated dibenzo-p-dioxins (PCDD) and places them in the dl-PCBs category 8 .Dioxin-like PCBs, known for their binding affinity to the aryl hydrocarbon receptor (AhR), may impair mitochondrial function via the AhR pathway and lead to diabetes 46 .Although the study did not directly investigate these mechanisms, the strong association with dioxin-like PCBs suggests a potential role for diabetes, which may be due to the mentioned mechanisms.

Strengths and limitations
To the best of our knowledge, this represents the first historical cohort study to explore the association between PCBs and type 2 diabetes in Iran.Prior investigations assessing the relationship between type 2 diabetes and environmental risks of PCBs encountered some limitations, such as conflicting results, the type of studied population (often not representative of the general population), a limited number of diabetes samples, and inadequate consideration of background or confounding variables.In this study, we tried to reduce limitations by considering the level of occupational exposure, diet, underlying diseases of the participants, and the type of population (via selecting a sample from the general population).The main strengths of the present study are the increase in the probability of obtaining cases due to the ten-year follow-up period, which is a relatively long period.Also, diabetes diagnoses were based on fasting blood sugar tests (FBS) rather than solely on self-declaration.By the way, the study faced some limitations.First, due to the coincidence of the study period with the COVID-19 pandemic, it was not possible to include all diabetes individuals.Second, because of the impossibility of accessing the participants' medical records, some data for certain variables such as various diseases, the amount of aquatic food consumption, vitamin D, or other supplements was obtained based on their self-declaration in response to retrospective questions.Other limitations of the present study were due to various factors such as unwillingness to participate, death, relocation, or lack of access to participants.Also, as we used several measures of exposure to PCBs and we conducted multiple comparisons in our study, some significant relationships may was not beyond the chance finding and the corrected estimate of p-values may be bigger than reported values.This may be a limitation for p-values that were not far less than 0.05.For p-values that were extremely small, this limitation is not of great importance.

Conclusions
The study support a significant correlation between exposure to total PCBs, dioxin-like and non-dioxin-like PCBs, with type 2 diabetes after adjusting for potential confounders.These findings emphasis the importance and necessity of implementing educational, health, and policy activities to reduce PCBs exposure and prevent this disease.Additional research, including large longitudinal studies designed with samples from various cohort studies within a community (e.g., Iran), is recommended to thoroughly evaluate the effects and mechanisms by which influence type 2 diabetes, especially considering all the confounding factors of the disease.Such studies, with considering all potential confounders, enhance comprehensively the understanding of the findings of this study.Furthermore, the findings should motivate further research into the simultaneous examination of PCBs exposures and other reported environmental causal factors (e.g., Cadmium, PAHs, etc.) suspected of contributing to type 2 diabetes.Identifying and addressing the root causes and risk factors associated with the disease can lead to the development of more effective prevention and management strategies.Lastly, Bushehr is a coastal and industrial city, with a hot and humid climate, and a special food culture, so that fish and shrimp are the main food of the citizens; it is recommended to consider these issues in generalizing the findings of this study to other societies. Vol

Table 1 .
General characteristics of the study participants.HDL high-density lipoprotein, LDL low-density lipoprotein, SD standard deviation, BMI body mass index, FBS fasting blood sugar, IPAQ international physical activity questionnaire.a p < 0.2: for selecting covariates by simple logistic regression.

Table 2 .
Distribution of the PCBs serum concentration of the study population (n = 181).CI confidence interval, OR odds ratio, ΣPCBs total polychlorinated biphenyls, ΣNdl-PCBs non-dioxin-like-PCBs, Σdl-PCBs dioxin-like-PCBs, IQR interquartile range.a The numbers reported are the closest to the 50th percentile, and they were considered as an estimation of the median.b More than 70% of the samples were less than the detection limit, so the median value set half of the LOD.Vol.:(0123456789) Scientific Reports | (2024) 14:8562 | https://doi.org/10.1038/s41598-024-59308-9

Table 3 .
Multiple logistic regression analysis between type 2 diabetes and PCBs serum concentration.OR odds ratio, CI confidence interval, ΣPCBs total polychlorinated biphenyls, ΣNdl-PCBs non-dioxin-like-PCBs, Σdl-PCBs dioxin-like-PCBs.a The models were adjusted for age, waist circumference, family history of diabetes, and dyslipidemia.b p < 0.05 is significant.c p > 0.05: the model is fit.